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ABSTRACT 

Despite their unique astrophysical relevance, the outcome of white dwarf binary merg- 
ers has so far only been studied for a very restricted number of systems. Here we 
present the results of a survey with more than two hundred simulations systematically 
scanning the white dwarf binary parameter space. We consider white dwarf masses 
ranging from 0.2 to 1.2 Mq and account for their different chemical compositions. We 
find excellent agreement with the orbital evolution predicted by mass transfer stability 
analysis. Much of our effort in this paper is dedicated to determining which binary sys- 
tems are prone to a thermonuclear explosion just prior to merger or at surface contact. 
We find that a large fraction of He-accreting binary systems explode: all dynamically 
unstable systems with accretor masses below 1.1 Mq and donor masses above ~ 0.4 
Mq are found to trigger a helium detonation at surface contact. A substantial fraction 
of these systems could explode at earlier times via detonations induced by instabilities 
in the accretion stream, as we have demonstrated in our previous work. We do not find 
definitive evidence for an explosion prior to merger or at surface contact in any of the 
studied double carbon-oxygen systems. Although we cannot exclude their occurrence 
if some helium is present, the available parameter space for a successful detonation 
in a white dwarf binary of pure carbon-oxygen composition is small. We demonstrate 
that a wide variety of dynamically unstable systems are viable type la candidates. The 
next decade thus holds enormous promise for the study of these events, in particular 
with the advent of wide-field synoptic surveys allowing a detailed characterization of 
their explosive properties. 

Key words: supernova: general - white dwarfs - nuclear reactions, nucleosynthesis, 
abundances - hydrodynamics 



1 INTRODUCTION 

Type la super novae are lynchpins in measuring distances 
on cosmological scales (Riess et al. 1998; Perlmutter et al. 
1998). These events are bright enough to be seen from a 
great distance such that many thousands have been cata- 
logued by modern transient surveys such as PTF (Rau et al. 
2009), CfA (Hicken et al. 2009a,b), and LOSS (Leaman et al. 
2011). Despite the large number of detected events, the pos- 
sible progenitor (s) of type la supernovae are still not known 
(Howell 2011). 

A strong constraint comes from the type la event re- 
cently discovered in MlOl, which placed an upper limit of 
0.02 Rq on the size of star that exploded (Nugent et al. 2011; 
Bloom et al. 2011). This all but eliminates the possibility 
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that the companion was a red giant (Li et al. 2011), and 
restricts the pre-explosion mass transfer rate from any non- 
giant companion to a narrow range of values (Horesh et al. 
2011; Chomiuk et al. 2012). Observations indicate that the 
event is a fairly standard type la (Brown et al. 2011; Bloom 
et al. 2011), and thus the particular progenitor channel that 
led to this event is likely to be similar to a significant frac- 
tion of the type la events observed. The strong restrictions 
and/or elimination of non-degenerate companions compels 
us to explore the possibility that the companion is also a 
degenerate star. This double degenerate (DD) scenario for 
type la would also be consistent with the observed delay 
time distribution, which is approximately proportional to 
(Totani et al. 2008; Maoz et al. 2010; Ruiter et al. 2010). 
This progenitor channel would further explain why hydro- 
gen has never been observed in a type la (Howell 2011). 
Moreover, the recently observed super-Chandrasekhar mass 
events (Howell et al. 2006; Scalzo et al. 2010; Taubenberger 
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et al. 2011) give further credence to the idea that at least 
some fraction of type la may be attributed to the DD sce- 
nario. 

Because white dwarfs (WDs) represent the final evolu- 
tionary state of the vast majority of stars, they are also the 
most common constituents of binaries containing compact 
stellar remnants. WDs are repositories of thermonuclear fuel 
and it is not surprising that there are several pathways that 
can produce a powerful explosion. In fact, it seems that all 
but the lowest mass WDs (< 0.2 Mq) are capable of being 
ignited by impacts that are of order their own escape speeds 
(Rosswog et al. 2009; Raskin et al. 2009). Moreover, the 
resulting theoretical lightcurves and spectra show excellent 
agreement with observed type la. Unfortunately, these colli- 
sion events are most likely too rare to explain a substantial 
fraction of the observed type la events. 

More frequently, a WD will find itself as a companion 
to another WD after both stars have evolved off the main 
sequence or after a dynamical exchange in a dense stellar 
environment (Shara & Hurley 2002; Lee et al. 2010). Bi- 
nary population modeling, which is calibrated by the cur- 
rently observed systems, estimates around 10^ double WDs 
(DWDs) in the Milky Way (Nelemans et al. 2001b). The 
rates at which these systems come into contact, although of 
great interest for the formation of interacting binaries and 
possibly explosive mergers, remain uncertain, because they 
depend primarily on the distance of the two WDs after leav- 
ing the poorly understood common envelope phase (Woods 
et al. 2011). Once the WDs have exited this phase, the two 
WDs are driven together via gravitational wave radiation, 
the signature of which is expected to be a dominant com- 
ponent of the ambient gravitational wave sky (e.g. Lipunov 
et al. 1987; Hils et al. 1990; Nelemans et al. 2004; Liu et al. 

2010) . 

A critical stage in the evolution of a compact WD bi- 
nary is the period just after mass transfer has been initiated. 
For the binary to survive, the mass transfer must be stable, 
which depends sensitively on the internal structure of the 
donor star, the binary mass ratio, and the angular momen- 
tum transport mechanisms (e.g. Marsh et al. 2004; Gokhale 
et al. 2007). The torques applied by a disk that may form 
from accreted material, and the tides raised on the accretor, 
are crucial for the binary's orbital stability (Iben et al. 1998; 
Piro 2011). If these mechanisms are unable to widen the or- 
bit, the system will merge quickly. The subsequent evolution 
can produce a R Corona Borealis star (Webbink 1984; Clay- 
ton et al. 2007; Longland et al. 2011), trigger an explosion 
as a type la (Webbink 1984; Iben k Tutukov 1984; Piersanti 
et al. 2003a, b; Yoon et al. 2007), or undergo an accretion- 
induced collapse to a neutron star (Saio & Nomoto 1985; 
Nomoto k Kondo 1991; Saio k Nomoto 2004; Shen et al. 

2011) . 

While theoretical advances detailing the possible out- 
comes of WD mergers have been made in the last decade (see 
review by Postnov k Yungelson 2006) , simulations have only 
probed the parameter space of DD mergers for a very small 
number of systems, with the principle emphasis being on 
systems with a total mass greater than the Chandrasekhar 
mass. We have previously stressed the vital role of the nu- 
merical initial conditions in determining the outcome of a 




Figure 1. The initial chemical composition of double WD sys- 
tems with donor masses M2 from 0.2 to 1.05 Mq and accretor 
masses Mi between 0.2 and 1.2 Mq. See text for additional in- 
formation. 

merger (Dan et al. 2011), simulations with approximate ini- 
tial conditions start with too little angular momentum and 
predict a merger that occurs much more rapidly than what is 
realized in nature. Our previous results indicated that many 
DD systems with a donor composed largely of helium (He) 
are expected to be ignited by an interaction between the ac- 
cretion stream and the torus of material acquired from the 
donor (Guillochon et al. 2010). Here we greatly expand upon 
our initial study with 225 new simulations that fully explore 
the possible combinations of WD donors and accretors in the 
mass range between 0.2 and 1.2 Mq. We demonstrate that 
our systematic study exhibits the orbital stability trends ex- 
pected for merging DD systems, thus enabling us to evaluate 
which systems are expected to produce strong thermonu- 
clear events in the lead up to or at the point of merger. We 
find that many systems with He donors may produce explo- 
sions, but that systems with carbon-oxygen (CO) donors are 
incapable of exploding without the inclusion of a significant 
He layer. 

In Section 2, we recap our numerical approach and de- 
scribe how our simulations cover the parameter space of DD 
systems. In Section 3 we give a brief overview of how dynam- 
ical orbital stability is determined and compare these expec- 
tations to our simulation results. Section 4 discusses detona- 
tions that are either triggered by the accretion stream or at 
the moment of merger for both He and CO mass-transferring 
donor stars. We summarize our findings and relate our mod- 
els with type la observations in Section 5. For reference, we 
provide fitting formulae bounding the region of parameter 
space where detonations are expected in Appendix A. 



2 NUMERICAL APPROACH 

The simulations of this investigation are performed with a 
3D smoothed particle hydrodynamics (SPH) code (Rosswog 
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Figure 2. Remnant profiles both with accurate (sohd-red fine) and approximate (dashed-green fine) initial conditions (ICs) for the 
0.81 + 0.9 Mq. Shown are the spherically enclosed mass (left) and spherically averaged density (center) and temperature (right) vs. 
radius from the center of the accretor. 



et al. 2008), for recent reviews of the SPH method see e.g. 
Monaghan (2005) and Rosswog (2009). The orbital dynam- 
ics of a binary system is very sensitive to the redistribution 
of angular momentum during mass transfer and therefore ac- 
curate numerical conservation is key to a reliable simulation. 
SPH is a very powerful tool for this type of study because the 
equations conserve angular momentum by construction, see 
e.g. Section 2.4 in Rosswog (2009). In practice, the quality 
of conservation is determined by the force evaluation (in our 
case by the opening criterion of our binary tree Benz et al. 
1990) and by the time stepping accuracy, both of which can 
be made arbitrarily accurate. Our code uses an artificial vis- 
cosity scheme (Morris & Monaghan 1997) that reduces the 
dissipation terms to a very low level away from discontinu- 
ities, together with a switch (Balsara 1995) to suppress the 
spurious viscous forces in pure shear flows. The system of 
fluid equations is closed by the Helmholtz equation of state 
(Timmes & Swesty 2000). It accepts an externally calcu- 
lated nuclear composition and allows a convenient coupling 
to a nuclear reaction network. We use a minimal nuclear 
reaction network (Hix et al. 1998) to determine the evolu- 
tion of the nuclear composition and to include the energetic 
feedback onto the gas from the nuclear reactions. A set of 
only seven abundance groups greatly reduces the computa- 
tional burden, but still reproduces the energy generation of 
all burning stages from He burning to NSE accurately. We 
use a binary tree (Benz et al. 1990) to search for the neigh- 
bor particles and to calculate the gravitational forces. More 
details about our SPH code can be found in Rosswog et al. 
(2008) and the provided links to the literature. 



2.1 Initial conditions 

The WD donor mass distribution in our parameter space 
starts at 0.2 Mq and reaches up to 1.05 Mq in steps of 0.05 
Mq. For accreting WDs we investigate those with masses 
between 0.2 and 1.2 M©, again in steps of 0.05 Mq. Accu- 
rate initial conditions (ICs) were constructed following the 
procedure described in Dan et al. (2011), here we only briefly 
summarize it. 

For all our simulations the stars are initially synchro- 
nized, cold and isothermal (T = 10^ K). In a first step, the 
stars are individually driven to an accurate numerical equi- 



librium by adding a velocity-dependent damping term to the 
momentum equation. Subsequently, we place the stars in the 
corotating frame at an initial separation that is large enough 
to avoid any immediate mass transfer. The orbital separa- 
tion is then adiabatically reduced in a way that the timescale 
over which the separation changes is much longer than the 
dynamical timescale of the donor. During this process, spuri- 
ous oscillations are minimized by applying a damping force. 
When the first SPH particle reaches the inner Lagrangian 
point Li, we stop the relaxation process, transform the ve- 
locities to the fixed frame and set this moment as the time 
origin (t = 0) of the simulation. Most of the previous SPH 
WD merger calculations (Benz et al. 1990; Segretain et al. 
1997; Guerrero et al. 2004; Yoon et al. 2007; Loren-Aguilar 
et al. 2009; Pakmor et al. 2010, 2011) started from rather 
approximate ICs. In Dan et al. (2011) we have shown in 
detail that the initial conditions have a strong impact on 
all major aspects of the simulations. Here, we only briefly 
illustrate the differences between the approximate and our, 
accurate ICs in the remnant structures resulting from the 
merger of a 0.81 + 0.9 Mq system in Figure 2. When using 
the accurate IC there is less mass in the disk and more mass 
in the trailing arm than when starting from the approximate 
IC. The approximate IC also overestimates the density and 
temperature in the region surrounding the accretor. As we 
explain in Section 5, the pollution of the ambient environ- 
ment prior to any supernova-like event may have important 
observational consequences. 

The compositions used in this study are shown in Fig- 
ure 1. Below 0.45 Mq the stars are made of pure He. In 
principle, WDs with masses below 0.45 Mq could harbor 
heavier elements in their cores (Iben & Tutukov 1985; Han 
et al. 2000; Moroni & Straniero 2009), but Nelemans et al. 
(2001b) have shown that the probability for hybrid WDs in 
this range is 4 — 5 times smaller than for He WDs. Based on 
Rappaport et al. (2009), who found that a 0.475 Mq WD 
has a He envelope of about 30% of its total mass, between 
0.45 and 0.6 Mq we adopt a hybrid He-CO composition 
with a ~ 0.1 Mq He mantle and pure CO core. WDs be- 
tween 0.6 and 1.05 Mq are assumed to be made entirely of 
CO. According with Iben & Tutukov (1985), there is only a 
small amount of He at the surface of the CO core for this 
range of masses. It has been shown that tiny amounts of He 
can be important for detonations of CO matter (e.g. Seiten- 
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Figure 3. Number of orbits of mass transfer for the He (left) and CO (right) donors. For He mass-transferring systems, the mass transfer 
stabihty hmits are over-plotted. The region below the continuous lower line indicates where a disk forms, above this line the accretion 
stream directly impacts on the accretor. The mass transfer is guaranteed to be unstable above the upper continuous line (Marsh et al. 
2004). Dashed lines represent stability limits assuming a synchronization timescale of the accretor with the orbit of 0.1, 10 and 100 yr. 
All systems with CO donors (right) belong to the unstable, direct impact regime of mass transfer. 



zahl et al. 2009), but for the numerical resolution that we 
can afford in this large parameter study (20,000 SPH parti- 
cles/star) we cannot properly resolve the outer He layers. For 
this range of masses there should be more Oxygen than Car- 
bon (Iben & Tutukov 1985), and we have chosen mass frac- 
tions X{^'^C) = 0.4 and X{^^0) = 0.6 uniformly distributed 
through the star, similar to Loren-Aguilar et al. (2009). For 
WDs with masses above 1.05 M© we use X{^^0) = 0.60, 
X(^°Ne) = 0.35 and X(^^Mg) = 0.05 (similar to the com- 
position found by Gil- Pons & Garcia-Berro (2001) for a 1.1 
Mq WD), distributed uniformly throughout the entire star. 



3 ORBITAL STABILITY 

After one or more common envelope phases, a white dwarf 
binary can emerge close enough for gravitational radiation to 
drive the binary towards shorter orbital periods until mass 
transfer is initiated. This occurs at orbital periods between 
2 and 15 minutes depending on the characteristics of the 
donor. For the binary to survive, the mass transfer must 
be stable, which depends sensitively on the binary mass ra- 
tio and the size of the donor star. Dynamical instability is 
guaranteed for systems in which the mass ratio of the donor 
star to the accretor star q = M2/M1 exceeds 2/3 (see e.g. 
Nelemans et al. 2001a). For q < 2/3, the system can still sur- 
vive as an interacting binary if orbital angular momentum 
losses are balanced by other mechanisms. In many instances, 
the formation of a tidally truncated accretion disk can effi- 
ciently return the advected angular momentum to the orbit. 
The system can then survive and evolve to longer periods 
as a result of mass transfer. On the other hand, mass trans- 
fer can proceed directly into the surface of the accretor. As 
a result, the system is destabilized by transforming orbital 
angular momentum into the accretor's spin. A direct impact 
will occur if the radius of closest approach of the accretion 
stream rmin is less than the radius of the accretor Ri. An 
estimate of rmin is given by (Lubow & Shu 1975; Nelemans 



et al. 2001a) 

— = 0.04948 - 0.03815 log q 

a 

+0.0475 log^ q - 0.006973 log^ q, (1) 

where a is the orbital separation. The two very-short period 
binaries HM Cnc and V407 Vul belong to this so called direct 
impact category as they both satisfy rmin < Ri (Marsh & 
Steeghs 2002; Ramsay et al. 2002). 

Tidal coupling can also return angular momentum to 
the orbit, heating both WD members in the process. This 
tidal heating can influence the donor's mass transfer re- 
sponse and, as a result, affect the stability of the system. 
For the coupling to effectively stabilize the system, the syn- 
chronization torques must act on a timescale shorter than 
the timescale over which the rate of mass transfer exponenti- 
ates. Some WD binaries clearly avoid merger, as we see their 
survivors in the form of AM Canum Venaticorum (AM CVn) 
binaries (see reviews by Warner 1995; Nelemans 2005; Sol- 
heim 2010). Survival is however not guaranteed even when 
mass transfer is expected to be stable. Mass transfer can, for 
example, exceed the Eddington rate or spark a thermonu- 
clear explosion that might destroy the binary. 

In Figure 3 we show the number of orbits each system 
survives prior to merger. The mass transfer rates we are able 
to numerically resolve are so far above the Eddington limit 
that photons are effectively trapped and, as a result, no sig- 
nificant mass loss is anticipated to occur. While all systems 
in our survey merge, we expect that binaries belonging to 
the disk accretion regime would survive if we were able to 
resolve the initial rate of mass transfer. In our simulations, 
the degree to which the donor star overfills its Roche lobe 
is artificially enhanced by our finite resolution. This results 
in a mass transfer rate that increases too quickly and pro- 
duces an under-massive accretion disk that cannot return 
sufficient angular momentum to the orbit to avoid merger. 

Additionally, dynamically unstable systems are ex- 
pected to be desynchronized at the initial separations where 
numerically resolvable mass transfer sets in (New & Tohline 
1997). Constructing such initial conditions is challenging as 
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Figure 4. Temperature (left) and density (right) evolution over the last 8 orbital timescales for the 0.6 + 0.9 Mq system. Tmax is the 
maximum temperature, p(Tmax) its corresponding density and (T) is the SPH-smoothed temperature, see Eq. (2). (T) is computed using 
the particles within the smoothing length of particle with Tmax- Colored bands show zb 1, 2 cr (standard deviations) distributions around 
the mean values, the white line indicates the arithmetic temperature mean of the neighbors of the hottest particle. Time is measured in 
units of the initial orbital period, Torbit- 



desynchronicity is a strongly destabilizing effect. While all 
our simulations were started with synchronized stars, the ac- 
cretor becomes desynchronized after only a few orbital rev- 
olutions, and thus the systems quickly evolve to a state that 
is close to what is realized in nature. Despite the quick build- 
up to this desynchronized state, our simulations are able to 
capture the basic angular momentum transport mechanisms 
correctly as illustrated by the dramatic increase in the num- 
ber of orbits as we approach the disk accretion boundary. In 
comparison with previous work at higher resolution in which 
the white dwarf binaries merge within only a few orbits, our 
results illustrate the importance of carefully constructed ini- 
tial conditions. 



4 DETONATIONS AT CONTACT 

In this section we investigate which binary systems could 
produce immediate detonations at contact, first during the 
mass transfer, which in the following we refer to as "stream- 
induced detonations", and, if detonations from the stream 
interaction are avoided, during the final coalescence that we 
call "contact detonations". In Sec. 4.1 we briefly summarize 
the detonation criteria that we use in Sects. 4.2 and 4.3 
to discuss the detonation prospects for both He- and CO- 
transferring systems. 



4.1 Detonation criteria 

In its simplest form, a detonation is a shock wave that ad- 
vances supersonically into an unburnt, reactive medium. Be- 
hind the shock wave exothermic reactions take place that 
continuously drive the shock to wear down possible dissi- 
pative effects. Whether a detonation forms or not depends 
delicately on the exact local conditions. The detonation con- 
ditions that have been applied in previous work (Fink et al. 
2007; Ropke et al. 2007; Townsley et al. 2007; Jordan et al. 



2008) "all rather loosely decide whether detonation con- 
ditions are reached based solely on the peak temperature 
reached in the simulation above a certain density" (Seiten- 
zahl et al. 2009). Apart from local density and tempera- 
ture, the temperature profile and the size of the region be- 
ing heated determines whether a detonation is triggered or 
not. Depending on these factors, the critical radii have been 
found to vary between centimeters and hundreds of kilo- 
meters (Seitenzahl et al. 2009). In general, these values are 
substantially below the resolution lengths in the relevant re- 
gions that we can afford in this study (^ 10^ - 10^ cm) with 
its more than 200 simulation runs. 

Given the delicacy of the exact detonation conditions, 
we resort in the bulk of this work to a simple yet robust 
timescale argument. We assume that dynamical burning sets 
in when the thermonuclear timescale, Tnuc, drops below the 
local dynamical timescale Tdyn- When this happens matter 
heats up more rapidly than it can expand, cool and quench 
the burning. To compute Tnuc we use Tnuc = CpT/cnuc (e.g., 
Taam 1980; Nomoto 1982), where the nuclear energy gener- 
ation rate, Cnuc, is taken directly from the reaction network 
and Cp is the specific heat at constant pressure, which we 
calculate from the Helmholtz equation of state. For the lo- 
cal dynamical timescale we use Tdyn = H/cs, where Cs is 
the sound speed and H the local pressure scale height. We 
estimate the H as P/pg, where g is the local gravitational 
acceleration. 

In our SPH formulation, the energy variable that is 
evolved forward in time is the internal energy. The tempera- 
ture at each particle position is found at every time step by a 
Newton- Raphson iteration. This temperature Ta, a labelling 
the particle, is the best estimate for the local matter temper- 
ature, but small numerical fluctuations in the internal energy 
(say due to finite numerical resolution, finite SPH neighbour 
number, finite time stepping accuracy etc.) can lead to larger 
fiuctuations in the SPH particle temperature estimate. This 
is especially so in regions where matter is degenerate and es- 
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Figure 5. Particle temperatures within a cross-section slice with a thickness equal with four smoothing lengths for a 0.3 + 0.6 Mq 
system using 4 x 10^ (left) and 2 x 10^ (center) SPH particles. The right panel shows the probability distribution functions (PDF) of the 
temperature within the red circle (interaction radius of the particle with minimum {tuuc / T^yj^) (T)) shown in the left and central panels. 



sentially independent of temperature. To some extent, this 
could be alleviated by designing a Fehlberg-type integration 
scheme that rejects and re-takes time steps if their temper- 
ature changes are too large. In addition to being computa- 
tionally expensive, this approach may be ineffective in the 
case the fluctuations are not associated with the choice of 
time-step. Therefore, one should be very cautious in inter- 
preting temperatures of very few (substantially less than the 
typical neighbour number) hot particles, which may also be 
due to small internal energy fluctuations. 

Although we stress that Ta is the best local temperature 
estimate, we also calculate the SPH-smoothed temperature 
value 

{T)a = Y.—nWi\ra-n\,ha), (2) 

which is a robust lower limit. Here, pb and mt are the den- 
sities and masses of neighbour particles 6, W is the cubic 
spline kernel and ha the smoothing length of particle a. In 
the limit of infinite resolution, both temperature estimates 
would coincide. Although it is plausible that a detonation 
would be triggered when {tuuc / Tdyn) (T) < 1, we only con- 
sider a detonation as inescapable consequence beyond rea- 
sonable doubt if the more conservative ratio {tuuc / Tdyn) {{T)) 
is also below unity. 

As a matter of illustration we show in Figure 4 the evo- 
lution of temperature and density of a 0.6 + 0.9 Mq system 
during the last stages of its merger. The left panel shows the 
maximum particle temperature (green) together with the 
SPH-smoothed value (red) and the arithmetic mean value 
within the interaction radius (= 2h) of the hottest particle 
(white). The shaded regions indicate the 1 and 2 a devi- 
ations. The corresponding quantities for the densities are 
plotted in the right panel. To illustrate the geometry of the 
potential ignition region we show in Figure 5 (left and cen- 
tral panel) the particle temperatures within a slice around 
the orbital plane. For a fair comparison we chose the slice 
thickness in each case to be 4 times the smoothing length 
of the particle with the minimum ratio of (rnuc/Tdyn)(7^)- 
The interaction region of this particle is indicated in both 
panels by the red circles. The probability distribution func- 
tions (PDFs) of the particles in these interaction regions 
are shown in the right panel. Despite the very low numeri- 



cal resolution that can be afforded in this parameter study, 
the overall dynamics is captured accurately, see for example 
the excellent agreement of the numerical results with the 
analytical estimates for the orbital evolution shown in Fig- 
ure 3. Nevertheless, at the current low resolution the shown 
interaction regions can contain particles from different mor- 
phological parts of the flow. In the shown example the most 
promising particle has neighbours in the cold (T < 5 x 10^ 
K) incoming accretion stream, in the hot interaction region 
('-^ 2 X 10^ K) and in the moderately hot inner disk regions 
(- 10^ K). As a consequence the PDFs can be multi-peaked. 
As resolution increases and the neighbour particles come 
from a smaller and smaller local volume, the PDFs will be- 
come narrower, single-peaked and T and (T) will become 
increasingly similar. 

In the remainder of the paper we explore the Mi — 
M2 plane to determine the regions in which the detonation 
conditions we have described are fulfilled. 



4.2 Helium mass transferring systems 

As discussed in Section 2.1, some of the donors we consider 
are composed of He in their outermost layers. For those 
systems with either He cores or thick He shells, the mass 
transferred to the accretor will be pure He until either the 
donor is disrupted, or until the He layer on the donor is 
exhausted. As donors with a significant fraction of He may 
exist for 0.2 ^ M2 ^ 0.6 M©, the mass ratio of systems 
with He donors can vary quite widely, and thus both dy- 
namically stable and unstable systems with He donors are 
expected to exist. Those dynamically stable systems that 
are close enough to one another to initiate mass transfer 
will appear as AM CVn systems, whereas the dynamically 
unstable systems, which experience an exponential growth 
in the accretion rate, are only likely to be observed prior 
to contact. As WDs are supported by degeneracy pressure, 
the mass-radius relationship of WDs only weakly depends 
on the WDs' composition, and thus the ratio of He to CO 
in WDs is largely unknown, aside from those WDs that are 
too massive to contain any significant He. 
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Figure 6. Torus masses of the He donor systems. If ever Tnuc ^ 
'Tdym we measure the masses when the timescales for the first 
time become equal, otherwise we take the mass when the ratio 
Tnuc/Tdyn IS minimal. The white contours show Tnuc/Tdyn when 
the ratio is at a minimum, with the thick contour showing where 
the ratio is unity. A comparison with Figure 8 shows that a large 
fraction of WDs with He donors undergo accretion stream trig- 
gered He-detonations prior to contact. The torus masses at deto- 
nation can be as large as ~ 0.1 M©, but also significantly smaller 
for larger Mtot- Systems with a massive He donor and small disk 
masses at detonation may be capable of producing multiple sub- 
luminous events prior to final coalescence. 



4-2.1 He detonations in stable systems 

Observed AM CVn systems are primarily expected to ac- 
crete through a disk (Nelemans 2005), although a few sys- 
tems are likely undergoing direct impact accretion (Marsh 
k Steeghs 2002; Ramsay et al. 2002). While the tempera- 
ture necessary for dynamical He burning is somewhat higher 
than the virial temperature of a ^ 0.6 Mq WD, slow He 
burning can proceed in the accretor's acquired He shell over 
a timescale shorter than the thermal timescale, raising the 
entropy of the accreted material before the excess heat can 
be radiated away into space. This is in fact the reason why 
no single WDs with a significant He layer are expected for 
M > 0.6 Mq, as even if the WD begins its life with a layer of 
He, this layer will be converted into CO on a timescale that 
is significantly shorter than the Hubble time. If the density 
at the base of the He layer is large enough (i.e., if the He 
shell is massive enough) , the temperature eventually rises so 
much that the burning timescale becomes shorter than the 
convective timescale in the envelope, leading to a run-away 
process that produces a detonation (Bildsten et al. 2007). 

However, a detonation can only occur in these kinds of 
systems if the donor is massive enough such that it can pro- 
vide the minimum required mass. It turns out that the min- 
imum mass necessary is larger than ^0.2 Mq for accretors 
with mass Mi < 0.6 M© (Shen & Bildsten 2009), systems 
that may or may not be dynamically stable depending on 
the synchronization timescale. 



4-2.2 Stream-induced detonations prior to contact 

As we described in Guillochon et al. (2010), detonations of 
the He torus can be triggered by the accretion stream sev- 
eral orbits prior to the merger event. Using our significantly 
larger set of accreting systems, we have provided an updated 
version of the figure originally presented in Dan et al. (2011) 
showing the systems that are expected to have surface det- 
onations triggered by the accretion stream, see Figure 6. A 
thermonuclear runaway will occur prior to the merging event 
if Tnuc/Tdyn ^ 1, a couditiou that depends on the thermody- 
namic state of the He torus acquired through the accretion 
process. The conditions for burning are more easily met for 
systems in which the temperature and density of the mate- 
rial acquired pre-merger is large to begin with, leading to 
conditions favorable for uncontrolled burning once the ma- 
terial is compressed by the dynamical interaction, whether 
by the stream or by the donor itself. 

The pre-interaction conditions in the He torus depend 
on three primary factors: the mass of the accretor, the com- 
ponent of the stream velocity normal to the accretor's sur- 
face at impact, and the evolution of M with time. For in- 
creasing accretor mass, the larger gravitational pressure at 
the accretor's surface leads to larger hydrostatic pressures, 
and thus greater equilibrium densities and temperatures. 
Conversely, an increasing accretor mass shrinks the accre- 
tor's radius relative to the circularization radius, resulting 
in more "disk-like" accretion and thus a less-concentrated 
torus that is more supported by rotation than fluid pressure. 
A shrinking accretor also reduces the normal component of 
the stream velocity relative to the accretor's surface, and as 
a result the degree of compression achieved as a result of 
the stream's impact onto the accretor's surface is more lim- 
ited. These effects are manifest as the larg 

Gr Tnuc /Tdyn ratios 
toward the highest accretor masses in Figure 6. 

A necessary condition for successful stream detonations 
is that the stream's ram pressure is considerably larger than 
the hydrostatic pressure of the torus, which is only possible 
if the mass of the stream is comparable in mass to the torus 
that has already accumulated. This implies that the sec- 
ond derivative of the amount of mass transferred, M, must 
be comparable to M itself. Therefore, stream detonations 
are only possible when an exponentially growing M can be 
sustained for several orbits, a condition that is typical of sys- 
tems with abrupt orbital evolution. Because of the stochastic 
nature of the triggering event, a detonation can be realized 
at any point in time where the He torus mass is larger than 
the minimum shell mass required for a successful detona- 
tion (Bildsten et al. 2007; Shen & Bildsten 2009; Woosley 
& Kasen 2011), which we find can precede the final merger 
event by many orbits. 



4.2.3 Contact detonations 

If detonations from the stream interaction are avoided, the 
system has another chance to initiate a He detonation dur- 
ing the final coalescence of the two systems. Except for WDs 
that are almost equal in mass, the dynamical interaction 
during the final coalescence involves the donor being ripped 
apart by tides, while the accretor remains relatively unaf- 
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Figure 7. Column density snapshots corresponding to the moment when the system first crosses (Tnuc/Tdyn)(^) < 1? or, if that does 
not occur, when (Tnuc/Tdyn)(r) is minimum. The ten snapshots shown here are representative of all DD chemical composition and mass 
combinations shown in Figure 1. In this parameter space we fully explore the possible combinations of WD donors and accretors and 
investigate their orbital stability and whether detonations prior or at the merger moment are possible. Systems with a mass ratio close to 
one (diagonal dashed-line) show a quick disruption, within about 20 orbits, while those away from this line show a slow depletion of the 
donor (ie. low M) over dozens of orbits of mass transfer and the formation of an extended atmosphere surrounding the accretor. We find 
that a large fraction of the systems with He donors (below the horizontal dashed-line) are expected to explode prior to or at the point 
of merger. In contrast, the CO accreting systems (above the horizontal dashed-line) are unlikely to explode at or prior to the merger. 



fected. Approximately half of the donor's remaining mass 
falls on top of the previously accreted He, but at a speed 
that's only a fraction of the escape speed from the accretor 
as the shredded donor is partly held aloft by its angular mo- 
mentum (Figure 7). Thus the final coalescence can be viewed 
as a low-speed collision (Rosswog et al. 2009). In simulations 
of WD-WD head-on collisions, collisions between pairs of 
WDs with mass as low as 0.4 Mq He WDs were seen to pro- 
duce detonations. The mutual escape speed of two 0.4 Mq at 



the point of collision is ^2G(Mi + M2)/{Ri + R2) = 3x10^ 
cm s~^ , which is approximately one-third of the escape speed 
from a 0.8 Mq accretor. 

At the point of disruption, half of the donor star falls 
onto the accretor as the donor's Roche lobe rapidly shrinks. 
Once this has occurred, the donor star can be viewed as ma- 
terial whose orbital dynamics are primarily determined by 
the surviving accretor. This means that the common cen- 
ter of motion, the barycenter, shifts from being located at a 
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Figure 8. Logarithmic ratio of the thermonuclear timescale Tnuc to the local dynamical timescale Tdyn at the moment of merger over a 
grid spanning the entire parameter space of He (left) and CO donors (right) masses. The ratio was computed using both the individual 
particle temperature T (upper panels) and the SPH-smoothed temperature (T) (lower panels, see eq. 2). Thermonuclear runaway is 
expected for configurations above the contours with Tnuc/Tdyn = 1- 



point lying between the two WDs to lying within the center 
of the accretor. From the plunging donor material's point of 
view, which originally moved on a circular orbit, the point 
around which it orbits moves, placing the material on a 
highly eccentric orbit with a large r. This results in a col- 
lision that is more violent, with a plunging speed that can 
be almost as large as the accretor's escape speed. In sys- 
tems with a large q, the bary center is usually already close 
to the accretor's center of mass, and thus there is no rapid 
acceleration within the corotating frame, and thus the final 
coalescence is much more gentle. This behavior is clearly 
visible in the right panel of Figure 9, which shows that the 
temperature at minimum {tuuc / Tdyn) (T) tends to increase 
with increasing q for a fixed total mass. 

The optimal conditions for ignition are reached when 
the temperature is at its peak, close to the time of the 
donor's disruption, provided that the accretor is below 1.1 
Mq and the donor is above 0.4 M©, see the upper-left panel 
of Figure 8. For all systems with Mtot < O.OM©, individual 
particle temperatures T of 10^ K and greater are reached in 
the torus formed via the mass transfer, but the density is 
too low to lead to a dynamical burning event. If we consider 
the more conservative smoothed temperature criteria, only 
systems with Mtot > I.SMq are capable of detonating at 
contact. Under both criteria, systems with ^ ^ 1 are more 



prone to produce detonations. Systems that do not satisfy 
the criteria for detonating either by the stream mechanism 
or at the point of merger may be progenitors of He-novae 
instead (e.g., Yoon & Langer 2004). 

The ratio Tnuc/Tdyn is very large in the right- lower cor- 
ner of Figure 8, left panels, close to the region of disk accre- 
tion. As described in Section 3, we expect that the systems 
should be dynamically stable, as they do not show any sign 
of a merger after 70 — 90 orbits of mass transfer. 

4.3 Carbon-Oxygen mass transfer 

While we find that He accreting WDs yield detonations, 
systems accreting CO are unlikely to explode at or prior to 
merger, see the right panels of Figure 8. Figure 9 shows the 
region of p — T plane where a CO detonation can be suc- 
cessfully ignited (Seitenzahl et al. 2009), together with the 
temperature, both T and (T), and density of all CO mass- 
transferring systems at the moment when {tuuc / Tdyn) (T) is 
minimal. Overall, there is a clear tendency to increase the 
temperature at minimum {tuuc / Tdyn) (T) with the mass ratio 
for a fixed total mass, see the right-panel of Figure 9. 

Based on the detonation criteria of Seitenzahl et al. 
(2009), Pakmor et al. (2011) have found detonations in their 
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logio P [g cm ^] Mtot 

Figure 9. Left panel: The filled diamonds show (p,T) at minimum {tuuc / T^yj^) (T) and are connected by black lines to the (p, (T)), 
represented by filled squares. Color coded is the total mass of the system in solar units. (T) is computed at the location of particle with 
temperature T. None of the CO mass-transferring systems reaches the CO detonation criteria of Seitenzahl et al. (2009) (shaded region) 
or the Tnuc/Tdyn < 1 Condition (see the right panels of Figure 8). Right panel: Temperature at minimum (tuuc / r^yn) (T) in q — Mtot 
plane. The temperature shows a clear trend towards higher values with increasing mass ratio and/or fixed total mass. 



merger study of (nearly) equal mass WDs. Because Pakmor 
et al. (2011) used approximate ICs that result in higher 
temperatures and densities during the evolution of mass 
transfer (Dan et al. 2011), we run a direct comparison with 
their 0.81 + 0.9 Mq system. We compared the minimum ob- 
tained value of Tnuc/^dyn with accuratc (Dan et al. 2011) 
and approximate (same as in Pakmor et al. 2011) ICs for 
non-rotating stars and found a difference of only 5%, with 
Tnuc/Tdyn ~ 7. Thc difference is much larger however, for 
systems with lower mass ratios. For the corotating system 
with 0.6 + 0.9 Mq components, Tnuc/rdyn for the approxi- 
mate ICs is four orders of magnitude greater than the value 
obtained for using the accurate ICs. The two ICs also lead 
to important differences in the merger remnant profiles at 
the moment when Tnuc/Tdyn is minimum: the temperature is 
overestimated when using the approximate ICs and the den- 
sity is larger in the hot envelope surrounding the accretor, 
where the dynamical burning takes place (Figure 2). 



5 DISCUSSION 

Although small patches of the white dwarf merger parame- 
ter space have been explored previously in a moderate num- 
ber of simulations (Benz et al. 1990; Rasio & Shapiro 1995; 
Segretain et al. 1997; Guerrero et al. 2004; D'Souza et al. 
2006; Yoon et al. 2007; Motl et al. 2007; Loren-Aguilar et al. 
2009; Pakmor et al. 2010; Fryer et al. 2010; Pakmor et al. 
2011; Dan et al. 2011; Raskin et al. 2011), the basic ques- 
tion: Which white dwarf systems produce explosions prior to 
or directly at merger? has so far remained unanswered. This 
serious gap in our understanding has in particular hampered 
our ability to confidently judge the potential of DD mergers 
as type la progenitor systems. In this paper we close this 
gap by methodically covering the whole parameter space of 
white dwarf mergers using more than 200 simulation runs. 
We have systematically scanned the white dwarf mass range 



from 0.2 to 1.2 Mq and accounted for their different initial 
chemical compositions. 

While such a broad study is necessarily limited in the 
numerical resolution that can be afforded in each simulation, 
we have found excellent qualitative agreement with previous 
numerical studies (Dan et al. 2011) and semi-analytical work 
(Marsh et al. 2004) on the orbital stability of white dwarf 
binary systems. One of the main conclusions of our previ- 
ous study (Dan et al. 2011) was that the accuracy of the 
numerical initial conditions has a major impact on all sub- 
sequent evolutionary stages and can be more important than 
numerical resolution. All our simulations were started from 
carefully constructed, initially tidally-locked binary configu- 
rations. The dynamically unstable systems in nature should 
-at the initial separations where numerically resolvable mass 
transfer sets in- have an accretor star that is desynchro- 
nized with the orbit (New & Tohline 1997). In practice, 
these initial conditions require a cumbersome approach as 
desynchronicity is a strong de- stabilizing effect, and there- 
fore constructing a relaxation scheme that does not lead to 
a premature merger is difficult. However, in all of our sim- 
ulations the accretor star becomes desychronized with the 
orbital period after about 5 to 10 orbital revolutions. We find 
systems that are near the dividing line between disk and di- 
rect impact accretion can still persist for many dozens of 
orbits prior to merging, despite the fact that the accretor is 
already desynchronous. We also find that as the binary ap- 
proaches merger, the orbital period changes so quickly that 
corotation cannot be maintained. In other words, while our 
initial conditions do not represent the physical state at the 
onset of resolvable mass transfer, the systems quickly evolve 
to the state that is likely realized in nature. 

As expected, systems of a given mass ratio lead to 
higher peak temperatures as the total system mass increases. 
We additionally find that there is a clear trend to produce 
higher peak temperatures as the mass ratio approaches unity 
(Figure 9). It comes as an interesting and maybe somewhat 
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Figure 10. Mass (in solar units) above a radius of 10-*^^ cm with respect to the accretor's center of mass for contact- (left) and stream- 
induced (right) He detonations. As in Figures 6 and 8, we calculate the masses when the timescales for the first time become equal, or, if 
that does not occur, we take the mass when the ratio Tnuc/Tdyn is minimal. The white contours show where the ratio Tnuc/Tdyn is unity. 



surprising result that a large fraction of the He-accreting 
systems are expected to explode prior to merger or at sur- 
face contact: we find that all systems with accretor masses 
below 1.1 Mq undergo such explosions, provided that the 
He-donating white dwarf exceeds ~ 0.4 (Figure 8). This 
opens up a large fraction of unexplored parameter space 
for thermonuclear explosions. A substantial fraction of these 
systems could even explode earlier via stream-induced deto- 
nations, as we have demonstrated in our earlier work (Guil- 
lochon et al. 2010). DD binaries made entirely of carbon 
and oxygen, in contrast, are unlikely to explode prior to or 
during the merger. This conclusion may need to be altered 
for cases where the CO white dwarfs possess more than a 
critical amount of He in their outer layers (Shen & Bildsten 
2009; Raskin et al. 2011). Systems that do not undergo an 
early thermonuclear event may still be prone to an explosion 
long after the merger (Yoon & Langer 2004). 

Based on a previous calculation (Fryer et al. 2010), it 
has been argued that the environment surrounding explod- 
ing DD systems is polluted with tens of percent of a solar 
mass and would have produced a detectable optical excess 
during the early evolution of 201 Ife (Nugent et al. 2011). 
However, in none of our exploding systems do we find more 
than 1% of a solar mass exterior to 10^° cm, and in fact 
many of our systems have no measurable mass beyond that 
distance (Figure 10). Because we are starting with initially 
synchronized WD binaries, the total amount of angular mo- 
mentum in our simulations is a robust upper limit. Therefore 
we can conclude that the ambient environment will contain 
even less mass at these distances than what our current sim- 
ulations predict. 

We have shown that DD mergers provide a broad chan- 
nel for producing thermonuclear explosions. While our sim- 
ulations indicate that systems undergo detonations within 
the He layer, it remains an open question whether this pri- 
mary detonation can actually trigger the explosion of the 
CO core. In the event that the CO core fails to detonate. 



these events may still resemble sub-luminous type la (Shen 
et al. 2010). Future research has to investigate which of these 
thermonuclear explosions are consistent with current type la 
SN observations. 



ACKNOWLEDGMENTS 

We thank Philipp Podsiadlowski and Ken Shen for insight- 
ful comments and stimulating discussions. We acknowledge 
support from DFG grant RO 3399 (M.D. and S.R.), the 
David and Lucile Packard Foundation (J.G., S.R., E.R.-R.), 
NSF grant: AST-0847563 (E.R.-R.), and the NASA Earth 
and Space Science Fellowship (J.G.). 



REFERENCES 

Balsara D. S., 1995, Journal of Computational Physics, 121, 
357 

W., Press W. H., Bowers R. L., 



Weinberg N. N., Nelemans C, 



Benz W., Cameron A. G 

1990, ApJ, 348, 647 
Bildsten L., Shen K. J., 

2007, ApJL, 662, L95 
Bloom J. S. et al., 2011, ArXiv e-prints 
Brown P. J. et al., 2011, ArXiv e-prints 
Chomiuk L. et al., 2012, ArXiv e-prints 
Clayton G. C, Geballe T. R., Herwig F., Fryer C, Asplund 

M., 2007, ApJ, 662, 1220 
Dan M., Rosswog S., Guillochon J., Ramirez-Ruiz E., 2011, 
ApJ, 737, 89 

D'Souza M. C. R., Motl P. M., Tohline J. E., Frank J., 

2006, ApJ, 643, 381 
Fink M., Hillebrandt W., Ropke F. K., 2007, A&A, 476, 

1133 

Fryer C. L. et al., 2010, ApJ, 725, 296 

Gil-Pons P., Garcia-Berro E., 2001, A&A, 375, 87 

Gokhale v., Peng X. M., Frank J., 2007, ApJ, 655, 1010 



© 2011 RAS, MNRAS 000, 1-13 



12 M. Dan, S. Rosswog, J. Guillochon, E. Ramirez- Ruiz 



Guerrero J., Garcia-Berro E., Isern J., 2004, A&A, 413, 257 
Guillochon J., Dan M., Ramirez-Ruiz E., Rosswog S., 2010, 

ApJL, 709, L64 
Han Z., Tout C. A., Eggleton R R, 2000, MNRAS, 319, 

215 

Hicken M. et aL, 2009a, ApJ, 700, 331 

Hicken M., Wood-Vasey W. M., Blondin S., Challis R, Jha 

S., Kelly R L., Rest A., Kirshner R. R, 2009b, ApJ, 700, 

1097 

Hils D., Bender R L., Webbink R. P., 1990, ApJ, 360, 75 
Hix W. R., Khokhlov A. M., Wheeler J. C., Thielemann 

F.-K., 1998, ApJ, 503, 332 
Horesh A. et aL, 2011, ArXiv e-prints 
Howell D. A., 2011, Nature Communications, 2 
Howell D. A. et aL, 2006, Nature, 443, 308 
Iben, Jr. L, Tutukov A. V., 1984, ApJS, 54, 335 
Iben, Jr. L, Tutukov A. V., 1985, ApJS, 58, 661 
Iben, Jr. I., Tutukov A. V., Fedorova A. V., 1998, ApJ, 

503, 344 

Jordan, IV G. C., Fisher R. T., Townsley D. M., Calder 

A. C., Graziani C., Asida S., Lamb D. Q., Truran J. W., 

2008, ApJ, 681, 1448 
Leaman J., Li W., Chornock R., Filippenko A. V., 2011, 

MNRAS, 412, 1419 
Lee W. H., Ramirez- Ruiz E., van de Ven G., 2010, ApJ, 

720, 953 

Li W. et aL, 2011, ArXiv e-prints 

Lipunov V. M., Postnov K. A., Prokhorov M. E., 1987, 
A&A, 176, LI 

Liu J., Han Z., Zhang F., Zhang Y., 2010, ApJ, 719, 1546 
Longland R., Loren-Aguilar P., Jose J., Garcia-Berro E., 

Althaus L. G., Isern J., 2011, ApJL, 737, L34 
Loren-Aguilar P., Isern J., Garcia-Berro E., 2009, A&A, 

500, 1193 

Lubow S. H., Shu F. H., 1975, ApJ, 198, 383 
Maoz D., Sharon K., Gal- Yam A., 2010, ApJ, 722, 1879 
Marsh T. R., Nelemans G., Steeghs D., 2004, MNRAS, 350, 
113 

Marsh T. R., Steeghs D., 2002, MNRAS, 331, L7 
Monaghan J. J., 2005, Reports on Progress in Physics, 68, 
1703 

Moroni P. G. P., Straniero O., 2009, in American In- 
stitute of Physics Conference Series, Vol. 1111, Amer- 
ican Institute of Physics Conference Series, G. Giobbi, 
A. Tornambe, G. Raimondo, M. Limongi, L. A. Antonelli, 
N. Menci, & E. Brocato, ed., pp. 63-66 

Morris J. P., Monaghan J. J., 1997, Journal of Computa- 
tional Physics, 136, 41 

Motl P. M., Frank J., Tohline J. E., D'Souza M. C. R., 
2007, ApJ, 670, 1314 

Nelemans G., 2005, in Astronomical Society of the Pacific 
Conference Series, Vol. 330, The Astrophysics of Cata- 
clysmic Variables and Related Objects, J.-M. Hameury & 
J. -P. Lasota, ed., p. 27 

Nelemans G., Portegies Zwart S. F., Verbunt F., Yungelson 
L. R., 2001a, A&A, 368, 939 

Nelemans G., Yungelson L. R., Portegies Zwart S. F., 2004, 
MNRAS, 349, 181 

Nelemans G., Yungelson L. R., Portegies Zwart S. F., Ver- 
bunt F., 2001b, A&A, 365, 491 

New K. C. B., Tohline J. E., 1997, ApJ, 490, 311 

Nomoto K., 1982, ApJ, 257, 780 



Nomoto K., Kondo Y., 1991, ApJL, 367, L19 

Nugent P. E. et aL, 2011, ArXiv e-prints 

Pakmor R., Hachinger S., Ropke F. K., Hillebrandt W., 

2011, A&A, 528, A117+ 
Pakmor R., Kromer M., Ropke F. K., Sim S. A., Ruiter 

A. J., Hillebrandt W., 2010, Nature, 463, 61 
Perlmutter S. et aL, 1998, Nature, 391, 51 
Piersanti L., Gagliardi S., Iben, Jr. I., Tornambe A., 2003a, 

ApJ, 598, 1229 
Piersanti L., Gagliardi S., Iben, Jr. I., Tornambe A., 2003b, 

ApJ, 583, 885 
Piro A. L., 2011, ApJL, 740, L53 

Postnov K. A., Yungelson L. R., 2006, Living Reviews in 

Relativity, 9, 6 
Ramsay G., Hakala P., Cropper M., 2002, MNRAS, 332, 

L7 

Rappaport S., Podsiadlowski P., Horev I., 2009, ApJ, 698, 
666 

Rasio F. A., Shapiro S. L., 1995, ApJ, 438, 887 

Raskin C, Scannapieco E., Fryer C, Rockefeller G., 

Timmes F. X., 2011, ArXiv e-prints 
Raskin C, Timmes F. X., Scannapieco E., Diehl S., Fryer 

C, 2009, MNRAS, 399, L156 

Rau A. et aL, 2009, PASP, 121, 1334 

Riess A. G. et aL, 1998, AJ, 116, 1009 

Ropke F. K., Woosley S. E., Hillebrandt W., 2007, ApJ, 

660, 1344 
Rosswog S., 2009, NewAR, 53, 78 

Rosswog S., Kasen D., Guillochon J., Ramirez- Ruiz E., 

2009, ApJL, 705, L128 
Rosswog S., Ramirez-Ruiz E., Hix W. R., Dan M., 2008, 

Computer Physics Communications, 179, 184 
Ruiter A. J., Belczynski K., Benacquista M., Larson S. L., 

Williams G., 2010, ApJ, 717, 1006 
Saio H., Nomoto K., 1985, A&A, 150, L21 
Saio H., Nomoto K., 2004, ApJ, 615, 444 
Scalzo R. A. et aL, 2010, ApJ, 713, 1073 
Segretain L., Chabrier G., Mochkovitch R., 1997, ApJ, 481, 

355 

Seitenzahl I. R., Meakin C. A., Townsley D. M., Lamb 

D. Q., Truran J. W., 2009, ApJ, 696, 515 
Shara M. M., Hurley J. R., 2002, ApJ, 571, 830 
Shen K. J., Bildsten L., 2009, APJ, 699, 1365 

Shen K. J., Bildsten L., Kasen D., Quataert E., 2011, ArXiv 
e-prints 

Shen K. J., Kasen D., Weinberg N. N., Bildsten L., Scan- 
napieco E., 2010, ApJ, 715, 767 

Solheim J.-E., 2010, PASP, 122, 1133 

Taam R. E., 1980, ApJ, 237, 142 

Taubenberger S. et aL, 2011, MNRAS, 412, 2735 

Timmes F. X., Swesty F. D., 2000, ApJS, 126, 501 

Totani T., Morokuma T., Oda T., Doi M., Yasuda N., 2008, 
PASJ, 60, 1327 

Townsley D. M., Calder A. C, Asida S. M., Seitenzahl I. R., 
Peng F., Vladimirova N., Lamb D. Q., Truran J. W., 2007, 
ApJ, 668, 1118 

Warner B., 1995, Ap&SS, 225, 249 

Webbink R. F., 1984, ApJ, 277, 355 

Woods T. E., Ivanova N., van der Sluys M., Chaichenets 

S., 2011, ArXiv e-prints 
Woosley S. E., Kasen D., 2011, ApJ, 734, 38 



© 2011 RAS, MNRAS 000, 1-13 



How the merger of two WDs depends on the mass ratio 13 

Yoon S., Podsiadlowski P., Rosswog S., 2007, MNRAS, 380, 
933 

Yoon S.-C, Langer N., 2004, A&A, 419, 645 



APPENDIX A: POLYNOMIAL FITTING 
FUNCTIONS FOR HE DETONATIONS 

Approximate formulas obtained by polynomial fitting of the 
contours lines with Tnuc/Tdyn = 1 are presented for the He 
mass-transferring systems (see Section 4.2). The contour 
(Tnuc/Tdyn) (7") = 1 showu lu thc uppcr-lcft panel of Figure 8 
is approximated to better than 4% by (in solar units) 

M2 = 21.065Mi -60.751Mf + 63.515Mi 

- 28.569Mi + 5.1209, (Al) 

with Ml between 0.46 and 1.13 M©. 

The conservative estimate for contact He detonations 
(tuuc / Tdyn) {{T)) — 1 (lower-left panel of Figure 8) is ap- 
proximated to within 3% with 

M2 = 39.43IM1 - 127.35Mf + 152.17Mi 

- 80.163Mi + 16.329, (A2) 

with Ml between 0.65 and 1.13 M©. 

For stream detonations (see Figure 6) , the cubic polyno- 
mial used to approximate the data between 0.88 ^ Mi ^1.2 
M© is: 

M2 = -16.318Mf + 53.568M1 - 57.809Mi 

+ 20.888. (A3) 

With this function we are slightly overestimating the number 
of systems that may lead to a detonation for Mi ^ 0.93 as 
the curve is not single- valued in this region. Above this value 
Equation (A2) approximates the data to better than 2%. 
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